Mott states under the influence of fermion-boson conversion: 

Invasion of superfluidity 
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I study the influence of fermion-boson conversion near Feshbach resonances on Mott states of 
Cooper pairs and demonstrate possible invasion of superfluidity. The quantum dynamics of Fermi- 
Bose gases is studied using both an effective coupled (7(1)® (7(1) quantum rotor Hamiltonian and a 
coupled XXZ ® XXZ spin Hamiltonian. I also point out two distinct branches of collective modes in 
superfluid states, one of which involves anti-symmetric phase oscillations in fermionic and bosonic 
channels and is always gapped because of fermion-boson conversion. 



Recent activities on atomic systems near Feshbach res- 
onances have generated much interest. By varying the 
two-body scattering length near Feshbach resonances, 
several groups have successfully achieved fermionic su- 
perfluid in a strongly interacting regimei*2i£i4, Superfluid 
states near Feshbach resonances turn out to be related to 
the BCS-BEC crossover studied long time ago 5 i 6 i 7 i 8 . This 
was recently pointed out by a few g r0U ps2ii2iii4i2iiii4. 
These groups have also made various efforts to incorpo- 
rate the two-body resonance between Cooper pairs (open 
channel) and molecules (close channel) explicitly in the 
many-body Hamiltonian and clearly illustrated an inter- 
esting role of atomic resonances. Near a Feshbach reso- 
nance the usual Cooper pairing amplitude and molecule 
condensate wave function are proportional to each other; 
the Feshbach resonance introduces an effective interac- 
tion between fermions, the interaction constant of which 
is 7j7 B /(2/x — v). Here jfb is the resonance width, v is 
the detuning energy of molecules and \l is the chemical 
potential of fermions. 

In this letter I am going to study a Mott state under 
the influence of Feshbach resonances, especially the ef- 
fect of fermion-boson conversion. A Mott state of bosons 
or of Cooper pairs appears whenever bosons or Cooper 
pairs in lattices are strongly repulsively interacting and 
if the corresponding filling factors are integers. One of 
important properties of a Mott state is its incompress- 
ibility, or a finite energy gap in its excitation spectra. 
And when hopping is renormalized to be zero due to 
repulsive interactions, the number of particles at each 
site is strictly quantized and discrete; particles are lo- 
cally conserved. I am going to show that in general Mott 
states are unstable with respect to fermion-boson con- 
version near Feshbach resonances. The primary reason 
is that fermions or bosons involved in resonating conver- 
sion are not conserved separately. So the conversion not 
only mediates an attractive interaction between fermions 
as realized before, but also, more importantly violates 
the local conservation law and introduces a novel mecha- 
nism to transport particles in some limit. The necessary 
condition for particles to be localized, or the number of 
fermions or bosons should be quantized might no longer 
exist. I find the invasion of superfluidity into a part of a 



parameter space where only Mott states are expected to 
be ground states if there are no Feshbach resonances. 

The model I employ to study this subject is an M- 
orbit Fermi-Bose Hubbard Model (FBHM). Consider the 
following general form of FBHM 



H — Hf + Hb + Hfb\ 

Hf = ~tf £ (fi v Ji^ + h.c.) 

</cZ>,7?,<7 



fe,77,£ fe 

H b = -t b ^2 ( b i b i + h -°-) 

<kl> 

+ - tyWkh + Vb/^nbkjnbk - 1); 
fe fe 

k . ') 

+Vbf^nbkhfk. (1) 

fc 

Here k, n and a label lattice sites, on-site orbits and 
spins; n = 1,2, ...M, a =T,j. f\ na (fk V a) is the 
creation(annihilation) operator of a fermion at site k, 
with on-site orbital energy e n and spin a. bt(pk) is 
the creation (annihilation) operator of a boson at site 
k (For simplicity, I assume there is only one bosonic 
orbital degree of freedom at each site). The fermion 
and boson number operators are, respectively, fibk — 
b\b k ,h fk = J2 v ,aftriafkri<T- tf and t b are hopping in- 
tegrals of fermions and bosons respectively and hopping 
occurs over neighboring sites labeled as < kl >. A is the 
attractive coupling constant in the Cooper channel which 
I assume to be much larger than the rest of couplings. 
Finally, Vf , Vb and Vfb are the strength of repulsive in- 
teractions between fermions, bosons, and fermions and 
bosons in the density-density channel^ One further as- 
sumes VbVf > Vf b to ensure that homogeneous states are 
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stable) . FBHMs similar to Eq. were previously applied 
to study Bose-Fermi mixtures in optical lattices^; most 
recently an FBHM with fermion-boson conversion was 
generalized to study the BCS-BEC crossover in lattices^. 
In the absence of the conversion term, FBHM consists of 
decoupled (attractive) Fermi-Hubbard model and Bose- 
Hubbard model; main properties of latter are knownlS. 

I first study the large-n limit where n b and ro/(<C M), 
the average numbers of fermions and bosons are both 
much bigger than unity. Because A is much larger than 
other coupling constants, the ground state of fermions 
for the on-site part of the Hamiltonian Hf should be a 
BCS state. This suggests that it should be convenient to 
work with the following coherent state representation, 



9o(n bk ) 



[exp(-^ bfc )6[]"« 



k ribk 



i Y[(v v + u v eK^{-i<j)f k )fl v ^fl vi )\vac 



> 



(2) 



Here v v , u v are the coherence factor in the BCS wave 
function which minizes the total on-site energy; go(nbk) 
is a unity for n max + n b > n bk > n b -n max , n max is much 
larger than one. These states form a low energy Hilbcrt 
subspace and are orthogonal in the limit which interests 
us, or < {<f>' fk }; {(/> bk }\{<f>fk}; {<t>fk} > is equal to zero if 
<j)fk 7^ 4>'f k or cf) bk ^ <j)' bk . At last, in the coherent-state 
representation one can show that fif k /2 = id/d<fif k , and 
h bk = id/d(f> bk ; or 



[-n/ fe ,exp(— i(f>f k >)] = 6 k>k > exp(-i0 /fc ), 
[n bk ,exp(-i(p bk >)] = S k:k > exp(-i0 bfc ). (3) 

In the subspace spanned by the coherent states defined 
above, I find that the effective Hamiltonian is 



H eff = -Jf E ^(^/fc " hi) + -± E^/ fc ~ n f) 

<kl> k 



-Jb ^ C0S (^f>fc _ M + V b ^(h bk - n b ) 2 

<kl> k 

- Tfs cos((j>f k - <$> b k) + Vf b (n bk - n b )(h fk - nff) 



The exchange couplings J/ , J b and Tfb can be estimated 
as 



Ur,V n ui,v'ri 



J b — n b t b , T FB = ^FBy/m ^ u n v n ; (5) 



E n = W (e l; — /i) 2 + Aq is the quasi-particle energy and 
Ao is the on-site BCS energy gap. Furthermore, nf and 
n b are functions of /x, v and Vf tb j b : 



2V b {ii + V f /A) - V bf (2fi -v + V b ) 



rib 



V f Vb 

f v b )- 



V 2 



V bf (2 f i + V f /2) 



2{VfV b 



V bf) 



(G) 



Obviously, the detuning energy v has to be sufficiently 
small in order for the population of bosons n b to be pos- 
itive. 

Eq.(j2|) and Eq.Q define the low energy quantum dy- 
namics of fermions and bosons under the influence of 
fermion-boson conversion near Feshbach resonances. In 
the absence of both Feshbach resonances (r fb = 0) and 
Vfb, the effective Hamiltonian describes two decoupled 
sets of quantum U(l) rotors in a lattice, the behaviors of 
which are well known. If n//2 or n b is a positive integer, 
the effective model can be used to study superfluid-Mott 
state transitions. A Mott phase corresponds to U(l) 
symmetry restored states and [/(l)-symmetry breaking 
solutions represent a superfluid phase. For the bosonic 
(Cooper pair) sector, the phase transition takes place 
when Tf — z.Jf/Vf {r b — zJ b /V b ) is equal to a critical 
value rf c (r bc ) (z(> 1) is the coordination number). The 
critical values which are of order of unity are usually cal- 
culated numerically. 

In the presence of Tfb, the Hamiltonian describes a 
coupled U(l) £§> U(l) quantum rotor model in a lattice. 
U(l) (g> U (1) symmetry breaking solutions when both r^ b 
are much larger than unity correspond to a superfluid 
phase. In the superfluid phase, the Hamiltonian in Eq.QJ 
further leads to the following semiclassical equation of 
motion in the long wave length limit 



d(j)f k 
dt 



dt 

1 dShf k 

2 dt 
dSn bk _ 

dt 



VfSh fk + 2V fb 5n bk , 

2VbSh bk + VbfShfk, 
= JfA(f>f k + T FB (4>bk - 4>fk), 
J b A<p bk + T FB {4>fk - 4>bk)- 



(J) 



Here 5hf k , b k = nfk.bk — Tij,b- I have taken a continuum 
limit and fc labels the coordinate of phases of bosons and 
fermion pairs (4>fk,bk) in this equation; A is a Lapla- 
cian operator. The lattice constant has been set to be 
one. This semiclassical equation shows that there are 
two branches of collective modes the dispersion relations 
of which are given below: 



a)uj 2 = a|Q| 2 , ^(w, Q -► 0) = Q 0); 
6V 2 -r»2+/3|Q| 2 , 

v f - v fb 



<j>f(w,Q -> 0) 



v b - Vfb 



(8) 



cj)f tb (u>,Q,) are the Fourier components of phase fields 
4>fk,bk(t)- In general, a and f3 depend on various 
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parameters in the Hamiltonian; and Qo is always pro- 
portional to TpB and ct on the other hand is indepen- 
dent of r FB . When V fb = 0, V> = V b = V and 
Jf = Jb = Jq, the dispersion relations are given as 



troduced two self-consistent order parameters 



2T FB V , a = (3 = 2J V . 



It is worth emphasizing that in the long wave length 
limit, mode a) is fully symmetric in phase oscillations of 
fermions and bosons, independent of various parameters; 
mode b) represents out-of-phase oscillations in fermionic 
and bosonic channel and becomes fully antisymmetry 
when Vb — Vf. In the absence of conversion (T FB = 0), 
these two modes correspond to two gapless Goldstone 
modes associated with breaking two decoupled U(l) sym- 
metries. However, in the presence of Feshbach resonances 
only the symmetric mode a) remains gapless correspond- 
ing to the usual Goldstone mode of superfluid while the 
antisymmetric mode b) is fully gapped because of the 
phase-locking effect of Feshbach resonances. 

In general, the wave functions for the many-body 
ground state and excitations ^ n ({(t>bk}', {<f>fk\){n = 
is the ground state) are the eigenstatcs of the Hamil- 
tonian in Eq.pjl, The boundary conditions are period- 
ical along the directions of 4>fk,bk with a period 2tt, so 
the wave functions are effectively defined on an S 1 ® S 1 
torus with radius of each S 1 equal to 2n. If the average 
number nf/2 and rib are integers, one introduces a gauge 
transformation ^ — ► \t Yik ex P( — ^ n f ( t > fk/2 — irib4>bk)\ the 
shifted number operators become 5hfk/2 = id/d(j)fk, 
Sfibk — id/dcf>bk- In this new basis, a spontaneous- 
symmetry breaking solution with wave function W ~ 
Uk^i^fk - 0o)<5(0bfc - </>o) represents a typical super- 
fluid state. A symmetry-unbroken solution with wave 
function ^ ~ fTfe ( 27r ) _1 exp(?m/ fe (/)/fe) ® exp(im b k4> m k) 
(jnfk,bk = for all k) on the other hand corresponds to 
a Mott state with 5ntk(hk)^ — or fifk(bk)^ — n f(b)^ a t 
each lattice site. 

To understand the influence of Feshbach resonances 
on Mott states, I now consider a situation where again 
both rib and nj/2 are integers and r/ is much less than 
r/ c so that Cooper pairs are in a Mott state in the ab- 
sence of Feshbach resonance. On the other hand r& is 
much less than the critical value r bc so that bosons are 
condensed. For simplicity, I have also assumed that Vbf 
is much smaller than Vt so that it can be treated as a 
perturbation. I am interested in the responses of Cooper- 
pair Mott states to fermion-boson conversion and carry 
out the rest of discussions in a mean field approximation 
(MFA). In this MFA, 4>jk = 4>f, 4>bk = <t>b for any lattice 
site k. The ground state ^fo(4>b, 4>f) (again defined on an 
S 1 <g) S* 1 torus with radius 2ir) is the lowest energy state 
of the following mean field Hamiltonian 



d 2 



i:) 2 



d d 



84>f d(f)b 
—z(JfAf cos0/ + JbA b cos</>b) 
-T FB cos((j) f - (f> b ). (9) 
Here again z is the coordination number; I have also in- 



A, 



-< cos <j)f >, Ab =< cos 4>b > 



(10) 



Here <> stands for an average taken in the ground 
state. Notice that the order parameters defined above 
are nonzero only when the U(l) symmetries are broken; 
particularly, Af is proportional to the usual BCS pair- 
ing amplitude. Following Eq.@ and discussions above 
one indeed shows that < ft v ^fl vl >= (52 v u v v v )Af , 

< bl >— y/nb~Ayl&. And A/ft vanish in a Mott state 
and are nonzero for superfluid states. 

As zJb is much larger than Vb, 4>b has very slow dy- 
namics; and the corresponding ground state for <pb can 
be approximated as a symmetry breaking solution. That 
is 



*o/(0/) = 



2tt' { Vf f 



I FB 

v f 



Aft) cos 0/111) 



Here At should be approximately equal to one in this 
limit; and in the zeroth order of Vf, \&o/ doesn't break 
the U(l) symmetry and stands for a Mott-state solution. 
Finally taking into account Ea.Q10[l and (|llfl . one finds 
that the self-consistent solution to A/ is 



; 2 V f 1 2 V f J 



(12) 



In the absence of T F b, Af vanishes as expected for a 
Mott state. However, the Mott state solution is unstable 
in the presence of any Feshbach conversion and the pair- 
ing order parameter A/ is always nonzero in this limit 20 . 

I want to emphasize that the average number of 
fermion per site is not affected by the fermion-boson con- 
version and remains to be an even integer (2/); rather, 
closely connected with the instability is the breakdown 
of particle quantization. Indeed, one obtains in the MFA 
the following results for hfk, < fifk >— rif = 21, 
< 5 2 h f k >= 1/2(T FB /Vf) 2 m 2A 2 . This illustrates that 
the resonance between states with different numbers of 
fermions at a lattice site eventually leads to a nonzero 
pairing amplitude A / . 

However, if rbj are much smaller than Tbcjc, then both 
Cooper pairs and bosons are in Mott states. Similar cal- 
culations lead to self-consistent solutions A/ = Aft = 

- \ i i / 1 iiimvm inmiM-l -mill' I 1 w i I i \ I 1 i w \ - 1 1 i i . - / - , \ i ■ i ■ / 1 1 - i i ■ l / \ i i . 1 l/ii 1 

= < COs(0ft ± (j)f) 



.19 



A" - Tfb A+ - 

2{v f + v b y bf ~"- 



(13) 



The second equality above simply shows the absence of 
superfluidity. But the first one indicates a subtle hidden 
order in the Mott states under consideration. Notice that 
A^* ~< bf.fk v ifkrjl > an d represents a tri-linear order. 
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The main conclusions arrived here do not depend on 
the large-n approximation introduced here. One can con- 
sider the opposite limit by assuming M = 1 and there 
is only one orbital degree of freedom at each lattice site. 
The two interaction terms (with two interaction constant 
A and Vf) in Hf (see Eq. lfl]l) in the single orbital limit 
can be rewritten in one term: V^/4^ fc hfkififu — 1) if 
one identifies Vj — Vf — 2A. In the limit where A is much 
larger than Vf, fermions are paired at each lattice site. 
Furthermore I assume bosons have hard core interactions 
(V b = oo) such that there can be only zero or one boson 
at each site. 

So the low energy Hilbert subspace Sk at each lattice 
site k consists of four states: 1) no Cooper pair, no bo- 
son; 2) no Cooper pair, one boson; 3) one Cooper pair, no 
boson; 4) one Cooper pair, one boson. They also corre- 
spond to a product of two pseudo spin 5=1/2 subspaces: 
St = S fk ® S bk , \<i z fk = ±1 >e Sf k , = ±1 >e S bk ; 
Sk is the on-site Hilbert space, and Sfk, b k are the on-site 
pseudo spin spaces for fermions and bosons respectively. 

Wfk = 1 >= fl]fli\ vac >/> Wbk = - 1 >= \ vac >/; 

Hk = 1 >= b l\ vac >b, Wlk = -1 >= \vac > b . 
\vac >f,b are the vaccuum of fermion and bosons re- 
spectively. Finally, in this truncated subspace, the fol- 
lowing identities hold a+ k = f^f^, aj k = fklfkt, 

a fk = 2 ffk-tffkjfklffk1 - 1; ct+ = b\, a^ = b k , 
al k = 2b\.bk — 1. So to have superfluidity, either a bk 
or (Jfk, or both of them need to have a finite expec- 
tation value in the XY plane. For instance to have 
fermionic superfluid, the expectation values of a^ k need 
to be nonzero. 

The effective Hamiltonian can then be written as 1 - 

H hs = - J l E + « - K E °h 

<kl> k 

- J } E °fk°fi + ° v fk° v f i - °fk°h - h ) E a h 

<kl> k 



~^FB E a fk°bk + a % a bk- (14) 
fe 

Here Jj = t)/V' f , J 1 = t b , T FB = lFB , h) = y. + V}/2 
and = /i — v/2. The Hamiltonian is invariant under 
a rotation around the z-axis or has an XY symmetry. 
The z-direction fully polarized phase of pseudo spins a b k 
((Jfk) represents the Mott phase of bosons (fermions) and 
the XY symmetry breaking states of pseudo spins a bk 
((Jfk) stand for the superfluid phase of bosons (fermions). 
The fermionic sector of this Hamiltonian was previously 
obtained and studied^i; it was also used to study BEC- 
BCS crossover in lattices^. 

When Tp B = 0, the Mott phase of Cooper pairs with 
filling factor equal to one occurs when is much larger 
than Jf . Assume that in this case is much less than J b 
so that a b k are ordered in the xy plane; then bosons are 
condensed. Taking into a finite amplitude of one 
considers a solution where the pseudo spin symmetry of 
a b k is spontaneously broken along a direction in the XY 
plane specified by < a bk > (the expectation value is taken 
in the ground state). In the molecular mean field approxi- 
mation, the effective external field acting on pseudo spins 
afk is h/, e // = zJf < a fk > +r^ B < a bk > +h z f e z ; and 
< afk > is calculated self-consistently in the ground state 
when hf. e ff is applied. One then arrives at the following 
solution < a fk > ■ < a b k > = ^fb/^ 1 } wnere < <Jfk > 
has been projected along the direction of < a b k > which 
lies in the XY plane. As mentioned before, development 
of such a component signifies superfluidity, or nonzero 
pairing amplitude. 

In conclusion, I have shown that certain Mott states 
are unstable with respect to the resonating fermion-boson 
conversion; in general superfluidity invades Mott phases 
because of the fermion-boson conversion. This work is 
supported by a Discovery grant from NSERC, Canada 
and a research grant from UBC. I acknowledge useful 
discussions with I. Affleck, E. Cornell and R. Hulct. 
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